Formation of Hydrogen Impurity States in Silicon and Insulators at Low Implantation 
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The formation of hydrogen-like muonium (Mu) has been studied as a function of implantation 
energy in intrinsic Si, thin films of condensed van der Waals gases (N2, Ne, Ar, Xe), fused and 
crystaUine quartz and sapphire. By varying the initial energy of positive muons (p"'') between 1 and 
30 keV the number of electron-hole pairs generated in the ionization track of the can be tuned 
between a few and several thousand. The results show the strong suppression of the formation of 
those Mu states that depend on the availability of excess electrons. This indicates, that the role of 
H-impurity states in determining electric properties of semiconductors and insulators depends on 
the way how atomic H is introduced into the material. 
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The implantation of energetic (MeV) positive muons 
(/i"^) in insulators or semiconductors commonly leads to 
the formation of the hydrogen-like bound state muonium 
[Mu = (/i+e^)] with a final charge state which can be ei- 
ther positive (Mu+), neutral (Mu°), or negative (Mu^). 
In semiconductors Mu is used to identify and investi- 
gate the electronic properties and the behavior of isolated 
hydrogen-like states and hydrogen- related impu- 

rities which are of fundamental and technological interest 
due to their influence on the electrical and optical prop- 
erties. Isolated H atoms in materials are difficult to de- 
tect by other spectroscopic means, which mostly require 
high H concentrations (see 0, @ and references therein 
for examples of vibrational spectroscopy studies of Si). 
In contrast, Mu - which behaves like a light H isotope 
(to^ ~ mp/9) - can be easily detected and characterized 
by the muon spin rotation (/iSR) technique due to its 
high sensitivity per spin. Muonium states, formed after 
implantation of energetic /i+ , remain isolated during the 
observational time window of the order of the /i"*" life time 
(2.2 fj,s). Therefore, a large amount of experimental infor- 
mation on the formation, structure and electrical activity 
of isolated H states in semiconductors has been obtained 
from /iSR, which has played a pioneering role in the iden- 
tification and characterization of hydrogen-like states. In 
Si, Ge and semiconductors of the III-V family two Mu 
states lying deep in the band gap have been identified at 
low temperatures (< 50 K)[l|: normal Mux in the tetra- 
hedral interstitial site with a large isotropic hyperfine in- 
teraction (hfi), and anomalous Mubc a-t a bond-center 
between two host atoms with a smaller, anisotropic hfi. 
In covalent semiconductors, Mut acts as an acceptor and 
Mubc as a donor. Recently, novel, very weakly bound 
Mu states {shallow Mu, binding energies 15-60 meV) 
with very low hfi have been established in a number of 
II- VI and III-V (nitrides) compounds [i,0,[l- Theoret- 



ical work has shown an universal alignment of hydrogen 
levels in semiconductors and insulators from which 
the electronic properties of hydrogen impurities can be 
derived. The predicted shallow donor hydrogen states 
in InN and ZnO have been confirmed experimentally by 
/iSR 

However, it has to be kept in mind that by the tech- 
niques used so far the spectroscopically investigated H- 
isotopes are energetically inserted in the solid. This re- 
sults in a large number N^h of electron-hole pairs gen- 
erated during slowing down of the incident particle. For 
instance, all /iSR experiments performed up to now used 
MeV-/i+ beams that generate 10^ - 10^ electron-hole 
pairs per implanted /i+ in the ionization track JOj. Sim- 
ilar or higher numbers of excess are created by the 
implantation of H or D ions used in the case of vibra- 
tional spectroscopy or channeling experiments. A sizable 
fraction of these electron-hole pairs escapes prompt re- 
combination and is still present around the thermalized 
impurity as shown in /iSR experiments with applied elec- 
tric field E of both polarities. The E-field clearly changes 
the Mu formation probabilit y by p ushing track and 
/i+ apart or together [lli[l2, 13, IJ], demonstrating that 
a significant fraction of Mu in semiconductors and insu- 
lators is generated by the capture of a track after the 
/i+ has stopped at an interstitial or bond site. In semi- 
conductors it appears that this so-called delayed Mu for- 
mation (in contrast to prompt Mu, where Mu forms dur- 
ing slowing down in charge-exchange cycles, followed by 
thermalization of Mu due to elastic collisions [31) is the 
origin of Mubc and for the recently discovered shallow 
Mu centers in III-V and II- VI semiconductors [l2l. [isl. ITst . 
The question therefore arises, whether and how the final 
states are influenced by the formation process, which is 
essential for studies on technologically important semi- 
conductors and insulators. This can be studied by using 
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the polarized low-energy /x+ (LE-zLt"*") beam at the Paul 
Scherrer Insitute (PSI, Villigen, Switzerland) [H, 17 1 
with variable implantation energy between 1 and 30 keV. 
It allows to investigate the formation of hydrogen-like Mu 
impurity states as a function of energy, i. e. as a function 
of Neh ■ By varying the energy, Neh can be tuned between 
a few and several thousand. This is up to five orders of 
magnitude less than for conventional MeV-muon beams. 
Below 1 keV nearly no track products are generated, thus 
approximating the case where H impurities are thermally 
introduced, which is the way, how trace atoms are incor- 
porated in the lattice in the course of wafer growth and 
fabrication processes. 

In this Letter we investigate for the first time the for- 
mation of thermal Mu as a prototype for isolated H im- 
purities as a function of implantation energy. In addition 
to intrinsic Si and sapphire (AI2O3) with more than one 
type of Mu we investigated thin films of van der Waals 
solids (s-Ne, s-Ar, s-Xe, S-N2) and fused and crystalline 
quartz (Si02) due to their simplicity concerning the fi- 
nal charge states: only one type of Mu*^ exists with an 
isotropic hfi close to vacuum Mu°. We find that delayed 
Mu formation is energy dependent in the keV range in all 
the investigated samples. Below ~ 10 keV the formation 
of those H impurity states that require a sizeable amount 
of excess e~ is strongly suppressed. The data on Si and 
AI2O3 support the interpretations that Mug^ in Si 
and Mu~ in AI2O3 [3| are formed by delayed capture of 
a track e~. 

The /iSR technique allows to differentiate between 
paramagnetic (Mu*^) and /i"*" in diamagnetic environment 
(free Mu+ or Mu~). Due to the hyperfine couphng 
between the /Lt+ and the e~ spin the observable Larmor 
precession frequency of isotropic Mu" is about 103 times 
larger than for the free /i"*". It splits into two intra-triplet 
lines that merge to one line at low fields (< 2 mT) where 
50% of the muon polarization is not observed due to un- 
resolved hyperfine oscillations between the triplet and 
singlet state. The diamagnetic and paramagnetic decay 
asymmetries and Amu were determined by measuring 
the amplitudes of the /x"*" and Mu precession signals in 
transverse (perpendicular to the /Lt+ spin) magnetic field, 
applied parallel to the sample normal. Ad and Amu are 
proportional to the fraction of muons in that particular 
states. 

The 0.5-mm thick Si sample with 50 mm diameter (un- 
doped, resistivity 10 kOcm, capped by a 2-nm thick oxide 
layer) was oriented with the (100) direction parallel to the 
sample normal. The quartz disc samples had thicknesses 
of 1 and 2 mm [Si02 crystal and fused quartz (Suprasil), 
respectively] and 50 mm diameter. The AI2O3 sample 
was a 0.5 mm-thick single crystal with 60 mm diame- 
ter. The solid gas films were grown at partial pressures 
between 10~^ and 5 x 10^^ hPa. Film thicknesses were 
about 1000 nm which is sufficient to stop all in the 

layer. At these deposition pressures grain sizes of order 
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FIG. 1: a) Typical jiSK asymmetry spectrum A(t) for the 
diamagnetic signal in s-Ar, and b) corresponding signal at 
low fields showing the 103-times faster Mu precession super- 
posed to the slow diamagnetic signal, c) Undoped Si, diamag- 
netic asymmetry Ad as a function of temperature T. Solid 
triangles: implantation energy Ei > 20 keV, open squares: 
El = 9.3 keV, open circles: Ei = 2.0 keV. d) Ad and Mut 
asymmetry a function of Ei. The lines in a) and b) 

are fits, and in c) and d) guides to the eye. 



100 nm are obtained ligl. For details on the experimental 
setup we refer to ref. [20|. 

Figure [T] shows typical /itSR asymmetry spectra, and 
displays the results for Si. In Si, at 5 mT only the pre- 
cession of the diamagnetic signal is observed. Due to 
limited statistics and time resolution Mu precession fre- 
quencies > 30 MHz are too high to be resolved with our 
present setup. In a field of 5 mT, the Mu^ intra-triplet 
lines are at about 70 MHz, whereas the Mugg transitions 
are between 35 and 50 MHz, depending on the orienta- 
tion of the B-field with respect to the (111) crystal axis. 
At 0.5 mT the MugQ frequencies are nearly unchanged 
and therefore not observable with our setup, whereas the 
7-MHz signal of Mu!p becomes visible. The 0.5-mT data 
are fitted with two components, a /i+ precession signal 
and the Mu^^ signal with exponential relaxation, whereas 
the 5-mT data are fitted with the /Lt"*" precession signal 
only. The temperature dependence oi Au at different im- 
plantation energies Ei is shown in Fig.[T^). Above 9 keV - 
corresponding to a mean implantation depth (d) of 65 nm 
and Neh — 2400 10] - Ad exhibits the same temperature 
behavior as in bulk /iSR experiments [l^ . The increase 
of Ad above ^ 150 K reflects the thermally induced ion- 
ization of MuJ^c- At = 2 kcV {{d) ^ 18 nm) and 
T < 150 K A/) is significantly larger than at higher ener- 
gies. As Fig.lTJj) shows the behavior oi Ad is not related 
to a change of the Mu^ fraction, which contrary to Ad 
does not depend on Ei. It rather reflects the unobserved 
ion which decreases with decreasing energy 
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FIG. 2: Muon and Mu asymmetries Ad and Amu as a function 
of implantation energy Ei for a) s-Ar grown at 6.5x10"'' hPa, 
and b) Si02 crystal and glass (Suprasil), T — 20 K. The 
magnetic fields are 10 mT (circles), 5 mT (triangles), 2 mT 
(stars) and 0.8 mT (squares). The lines are guides to the eye. 



and number of available track e" . This is also supported 
by the convergence of the two curves in Fig. [T]a,) at high 
T where Mug^ is ionized. The ~ 2-nm thick oxide layer 
present on the Si surface is too thin to explain the ob- 
served reduction of the Mu formation. Moreover, in a 
Si02 layer at low energy a higher Mu fraction should be 
observed, see Fig. [^b). The different dependence on the 
availability of excess e~ indicate that the main fraction 
of Mub(^ is due to delayed formation, whereas Mu!^ is a 
consequence of charge-exchange processes at epithermal 
energies - in agreement with bulk /iSR studies, where 
an applied E-field was used to vary the average distance 
between fi~^ and excess e~, and therefore the relative for- 
mation probability of these two states 

Figure [2] shows the energy dependence of Ad and 
^Mu for s-Ar (a) and Si02 (b). Only isotropic Mu is 
present, and and Amu represent a direct measure of 
the /i"*" and Mu fraction in the sample. The sum Atot 
= Ad + 2Amu ~ 0.263(1) is the total observable asym- 
metry, and there is within the experimental errors no 
missing fraction. Qualitatively, the s-Ar and Si02 data 
display the same behavior: with increasing energy Ad is 
decreasing while Amu is increasing correspondingly. The 
energy-dependent diamagnetic fractions Fd — Ad /Atot 
for various insulators are summarized in Fig.[3l With the 
exception of s-Ne all samples show a decreasing diamag- 
netic fraction with increasing energy. For Si02 and s-Xe 
bulk Mu fractions Fmu = (1 - Fd) of 85% and ~ 100%, 
respectively, are obtained at 20 keV [corresponding to 
(d) = 155 nm (SiOz), (d) = 185 nm (s-Xe)]. At this 
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FIG. 3; Comparison of the diamagnetic fraction Fd as a 
function of implantation energy Ei for different samples, 
B = 10 mT. Deposition pressures were; 7 x 10~® hPa for 
s-Ne, 6.5 X lO"** hPa for s-Ar, 2.2 x lO"'^ liPa for S-N2, and 
1.5 X 10~^ hPa for s-Xe. The lines are guides to the eye. 



energy the number of electron-hole pairs created in the 
ionization track is about 1000 [l3|. In the s-Ar and S-N2 
films even at the highest energy the observed Mu frac- 
tions (Fmu ~ 60%) are lower than the bulk results ob- 
tained with 4-MeV fi+ [Fmu - 100% (s-Ar), Fmu 80% 
(S-N2 at T < 30 K)]. The discrepancy is even more drastic 
for s-Ne where the film data are consistent with Fmu = 
in contrast to the bulk data with Fmu — 90% [14]. This 
disagreement can be explained by the suppression of Mu 
formation in granular s-Ne, s-Ar and S-N2 thin films, as 
we discuss below. 

The decrease of Fd with increasing Ei reflects the on- 
set of delayed Mu formation with increasing availability 
of excess e~ . From the flattening of Fd at ~ 20 keV 
we estimate the number of excess e~ necessary to satu- 
rate the delayed Mu yield to be of the order of thousand. 
The e~ may escape recombination with the /i+ by sev- 
eral processes: recombination with a cation from the ion 
track, trapping at grain boundaries, voids, and surfaces 
or escape from the surface (e~ escape depth 100 nm in 
s-Ar and s-Xe [llj, e~ mean free path in Si is ~ 20 nm 
at 300 K, increasing to > 100 nm at lower T). An ad- 
ditional obstacle for electron-muon recombination is also 
the large escape depth of 20 - 100 nm of epithermal /i"*" in 
wide band gap insulators such as S-N2 , s-Ar and s-Ne [l^] : 
after leaving the charge-exchange cycles where the last e~ 
are released the /z"'" may move such a distance away from 
its ionization track, further losing energy inefficiently by 
elastic collisions. This large e" — /i+ separation and the 
trapping of e" combine all together to suppress the de- 
layed Mu formation channel in s-Ne, s-Ar and S-N2. The 
total suppression of Mu formation in s-Ne is probably a 
consequence of a /i"*" escape depth larger than the typical 
grain size, making the formation of a delayed e~ — jX^ 
bound state unlikely. 

The energy dependence at T < 100 K of Ad in sap- 
phire (Fig. m shows an interesting anomaly compared to 
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FIG. 4: Diamagnetic asymmetry j4d as function of implan- 
tation energy Ei for sapphire, measured at B = 10 mT at 
different temperatures. The lines are guides to the eye. 



theoretical and experimental studies in this area. 

We are extending this kind of experiments to the in- 
vestigation of shallow Mu states: Preliminary data in 
ZnO show also a decreasing shallow Mu fraction with de- 
creasing implantation energy, which further supports the 
general trend found in other materials. 

This work was fully performed at the Swiss Muon 
Source S/xS, Paul Scherrer Institute, Villigen, Switzer- 
land. We are grateful to S.F.J. Cox for valuable discus- 
sions. We thank C. David from LMN at PSI for providing 
the Si sample. We thank M. Birke, Ch. Niedermayer and 
M. Pleines for their help in the initial phase of the exper- 
iment. The technical support by H.P. Weber is gratefully 
acknowledged. 



the data presented so far. At 100 K Ad decreases with 
increasing energy and reaches its smallest value of 0.025 
at 30 keV. This behavior correlates with the onset of de- 
layed formation of Mu" as seen in other insulators. The 
energy dependence oi Ad becomes less pronounced on 
reducing the temperature. At 4 K Ad exhibits a mini- 
mum at 10 keV and starts to increase again when further 
increasing the energy. This may reflect the delayed for- 
mation of diamagnetic Mu~, as suggested in a previous 
E-field fiSK experiment where the disappearence of Mu^ 
with increasing T is interpreted as thermal ioniziation of 
Mu^ with an activation temperature of 130 K [ll|. A re- 
cent theoretical work shows that could be the stable 
charge state in AI2O3 23]. Our data support this idea, 
and that Mu^ is formed by delayed e^ capture. 

In conclusion, the measured energy dependence of Mu 
formation in intrinsic Si and insulators shows as a gen- 
eral behavior that the formation of delayed Mu states re- 
quires the presence of the order of thousand excess e~ in 
the ionization track. With LE-/i+ H- impurity states can 
be studied without the generation of a non-equilibrium 
electron track. From the implantation energies involved 
we infer that the length scale of that part of the track 
that is involved in delayed Mu formation is of the order 
of 100 nm. At energies < 3 keV delayed Mu formation is 
nearly absent. This indicates that the formation of those 
H-impurity states which heavily depend on the availabil- 
ity of excess e~ is strongly suppressed in cases where the 
H-isotope is inserted in the solid without the concomitant 
presence of a sizeable number of excess e~ . This implies, 
that the role of H-impurity states in determining elec- 
tric properties of semiconductors and insulators depends 
on the way how atomic H is incorporated into the mate- 
rial. The question of the relative importance of different 
possible H-states and their occurrence as native impurity 
states in semiconductors and insulators is generally not 
addressed and we hope that our results will foster new 
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